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We study the problem of a magnetic polaron in an antiferromagnetic semiconductor (ferron). We 
obtain an analytical solution for the distortion produced in the antiferromagnetic structure due to 
the presence of a charge carrier bound to an impurity. The region in which the charge carrier is 
trapped is of the order of the lattice constant (small ferron) but the distortion of the magnetic 
structure extends over a much larger distance. It is shown that the presence of this distortion makes 
the ferron more stable, and introduces a new length scale in the problem. 

PACS numbers: 75.90.+W 

Keywords: Magnetic semiconductors, electronic phase separation, ferron, magnetic polaron, percolation 



I. INTRODUCTION. 

Quite a while ago Nagaev [J has shown that the min- 
imal energy for a charge carrier moving in an antifer- 
romagnetic background is obtained when the electron 
modifies the magnetic background and is self-trapped in 
a region with canted antiferromagnetic or ferromagnetic 
order. The term ferron (a magnetic polaron in an an- 
tiferromagnetic background) was coined there to name 
this new quasi-particle. Depending on the radius of the 
self trapping region, it is possible to differentiate between 
small ferron, localized in a region of the order of the lat- 
tice constant, and large ferron localized at larger scales. 
Small bound fcrrons are expected to be typical of the low 
doping region of the phase diagram of manganites, when 
the material is antiferromagnetic and insulating. 

There are also experimental data that confirm the exis- 
tence of bound ferrous in antiferromagnetic semiconduc- 
tors such as undcrcloped manganites. A review of them 
can be found in 0, Q . In reference Q , a liquid- like spa- 
tial distribution of magnetic droplets in Lai_a;Caa;Mn3 
with X — 0.05, 0.08 is reported. Ferromagnetic rich-hole 
droplets with a diameter of 4 — 5 lattice units isotropi- 
cally distributed with a mean distance of 9 lattice units 
among them are observed in an antiferromagnetic poor- 
hole background. Also it is reported that these mag- 
netic droplets are coupled together through the antifer- 
romagnetic background. In reference (3], a ^^Mn NMR 
study on the same compound within the doping range 
X < 0.23, confirms the electronic phase separation and 
report the existence of thin boundaries between antifer- 
romagnetic and ferromagnetic domains. These magnetic 
droplets contain a number of conduction electrons larger 
than one (about 30), but this is not essential in describing 
the spin distortion that they create in the antiferromag- 
netic background. 
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Almost in all the calculations of the fcrrons it was 
assumed that the region of the distortion of the mag- 
netic structure coincides with the range of localization of 
the electron. Often these regions were simply taken as 
spheres of radius R, to be determined self-consistently. 
However, as first pointed out by De Gennes 6], the dis- 
tortion of the magnetic order around a magnetic defect 
(i.e. ferron) may decay slowly with distance. In this 
paper we study what would be the "back effect" of this 
slowly decaying magnetic distortion on the conditions of 
the electron localization and on the properties of the re- 
sulting self-trapped state. First numerical calculations 
addressing this problem have been carried out recently 
by Nagaev ■ In this paper the magnetic anisotropy, al- 
ways present in real materials, was neglected, and he did 
not manage to obtain the results in a closed form and did 
not get proper estimates for the radius and the energy of 
the ferron taking this effect into account. 

In this article, we obtain the analytical solution for 
the distortion created by a localized conduction electron 
trapped in a region of the order of the lattice constant 
in an one-dimensional antiferromagnetic semiconductor. 
The main virtue of our calculation is than reveals the ex- 
istence of a region surrounding the trapping region which 
acts as an antiferromagnetic domain wall, as it was an- 
ticipated in reference ,7| . We calculated the properties of 
the ferron in this situation and show that the account of 
an extended magnetic distortion around it actually leads 
to an increase of the ferron stability and its binding en- 
ergy. Our result offers a key to understand the features of 
the phase diagram of antiferromagnetic semiconductors, 
such as manganites at the small doping range x < 0.1. 



II. RESULTS 

To study the problem of the magnetic polaron forma- 
tion, we start from the Hamiltonian proposed by Na- 
gaev 0. We analyze the one-dimensional case because 
it can be treated analytically and its solution provides 
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a full physical insight into the problem. Estimates to 
the three-dimensional case are given below. We consider 
an one-dimensional chain of magnetic ions along the y- 
axis, described by the double exchange (or Vonsovsky s- 
d) model. An uniaxial magnetic anisotropy term is added 
with X being the easy axis. The magnetic structure of the 
magnetic ions without the conduction electron is repre- 
sented by two sub-lattices with the spin up along the easy 
axis (that means = S) for ions at even sites {g = 2n) , 
and spin down (that means = —S) for ions at odd 
sites {g = 2n + 1). All distances are measured in terms 
of the lattice constant. A non-magnetic donor impurity 
is added to the chain at the point y = —1/2. Its con- 
duction electron is bound to it and it can jump between 
their two magnetic neighbors (bound ferron). This dis- 
turbs the pure antiferromagnetic order along the chain, 
even though the conduction electron is trapped only on 
these two magnetic ions. The Hamiltonian of the system 
can be represented by: 

Hsd = -t (^ati^„ao,a + a^ aa-i^a) 

-lY^S,S,+,-K'Y^{Slf (1) 

9 g 

where a^g., a^.o- are the conduction electron operator cor- 
responding to the site g and spin projection ct, sthe con- 
duction electron spin operator, Sg is the spin operator of 
the magnetic ion at site g (d-spins). The d-d exchange 
integral / is assumed negative in order to get the anti- 
ferromagnetic ordering. The anisotropy constant K' is 
considered positive. The energy of Coulomb interaction 
between the conduction electron and the impurity is an 
additive constant and, for this reason, omitted in the cal- 
culation. 

Depending on the relative value of the parameters W — 
2zt and AS, z being the number of first neighbors and 
S the magnitude of the d-spin, we have two different 
situations. In the case W » AS, we talk about a wide- 
band semiconductor. In the case W « AS, we talk 
about a double exchange semiconductor. 

Our goal is to obtain an expression for the magnetic 
energy of the system of d-spins, both in the case of wide- 
band and double exchange semiconductors. As 5 > 2 for 
the compounds of interest (typically, manganites), the 
d-spins are considered classically. Their orientations are 
described in a coordinate system centered in the position 
of each magnetic ion. As it was stated above, it is as- 
sumed that the conduction electron is in the lowest bound 
state in the space spanned by the operators of the sites 
g = —1,0. Also we assume that the magnetic moment 
of the ferron is directed along z-axis, that is, the con- 
duction electron spin acts as an effective magnetic field 
along the z-axis for the d-spin system. In this case, the 
following symmetries hold for the d-spin system: S^ = 0, 

and S^ = ^^(g+i). Then Sg = S {{-1)3 sm0g,O,cos0g). 



It is important to notice that here, although the d-spin 
lies in the z — x plane, the angle 9 is not the polar angle, 
but the spherical coordinate. 

To obtain the magnetic energy of the d-spin system 
the electronic part of the Hamiltonian must be aver- 
aged out. To do this, we assume that the electronic wave 
function is the ground state wave function for the domi- 
nant term in the electronic Hamiltonian, i.e. the hopping 
term in the wide-band case, and the s-d exchange term 
in the double exchange case. The other term is treated 
as a perturbation • In the case of wide-band semicon- 
ductor, the hopping term is diagonalised and its ground 

state is |<i>) — (0,-1 1/2 + ^01/2) 1^)- 

double exchange semiconductor, the ground state is the 

same but the operators being the operators with 

the spin projection along the direction of the vector Sg, 
instead of the laboratory z-axis. We treat the case of a 
wide-band semiconductor first. Then the magnetic en- 
ergy of the d-spin system is: 

jY,COS{9g+eg+l) 

a 

-L[cos0-i+coseo]- K^sin^Og-t (2) 

a 

where: J = -IS^, L = AS/ 4, and K = K'S^. 

Minimizing the equation Q with respect to the angles 
9g, a set of non-linear equations is obtained: 

J sin {Og + 9g+i) + J sin (6ig_i +9g)- 
-LsinOgiSg^^i + 5g^o] + K sm{2eg) = (3) 

There is a boundary condition 9g^±oo = Tr/2 if the 
chain is long enough, that means if KN » L, with N 
being the number of magnetic ions of the chain. Further, 
the above symmetry conditions imply that only the sites 
with g > must be considered. 

For the double exchange case, the same set of equations 
is obtained but the term in L in equation (0) must be 
changed by the standard effective hopping of the double 

exchange model, —tcos {^—^y^^ JJJj and the constant 
term —t must be changed by —2L. Using the symmetry 
condition 9-i — 9o, this is equivalent to change Lhy t/2 
in equation Therefore both cases can be treated on 
the same footing. 

Following the paper by Neel 8] , we look for a differen- 
tial equation for the d-spin distortion. If we were dealing 
with a ferromagnetic d-d exchange, the set of equations 
would describe a domain wall. To obtain a differen- 
tial equation for our problem an additional step is needed. 
Instead of working with the angles 9g, we perform a ro- 
tation of an angle tt around each a;-axis if the site is an 
odd site, and no rotation if the site is an even site. This 
corresponds to make the following changes in the angles: 
92n 92n, ^2n+i ^ TT - 02n+i- We assumc that the 
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length of variation of the angle 9 is larger than the lat- 
tice constant. We treat 9 as a, continuous function over 
the y-axis and perform a power expansion in the lattice 
constant. Then a differential equation is obtained. Tak- 
ing into account the above symmetries, we have to solve 
only for the positive semi-axis. Further, we divide the 
problem in two parts. For y > 0, we have: 



J--^+Ksm29^0 (4) 
For y — and using again the symmetries, we obtain: 



J — 

dy 



+ {J + K) sin 26*0 - L sin 6*0 = (5) 



where 6*0 = 9 {y = 0). This is the sine-Gordon equation 
with a boundary condition at the origin. 

We solve equation Q with the boundary condition 
that 9y-,+oo 7r/2. Multiplying by 2^ and integrating 
once, we obtain: 
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FIG. 1: Magnetization along z-axis, Sg, (in units of S) with 
L = 3 and anisotropy K = 2.5 ■ 10~^ (both in J units). Ana- 
lytical calculation. 



with a, b being two real constants to be determined. Note 
that a, 6 > guarantees that the angle 9 lies in the range 
(0,7r/2) for y belonging to the domain [0, -t-oo). Also 
note that at the infinite 9 goes to 7r/2, as required. The 
differential equation is identically satisfied if: 



J 



Kcos29 = C 



(6) 



where C is the integration constant. Following reference 
M, we made: 



and: 



8^ 
J 



(11) 



C 



K cos 29^ 







= J [Lsin6'o - (J + if)sin26'o]^ -/^cos26lo (7) 

The last step is to take into account equation ((Sj. Now 
we made the following change: 



smf 



/(y)-i 



/(y) + i 

The differential equation ^ is cast into the form: 



{f + lf f\dy) + J ' 



(8) 



L fo-l ^ J + K fo-1 
J fo + ^ J (/o + 1)' 



H 2 



//o-l 



J V./0 + 1 



(9) 



where we define fo ^ f [y = 0). The solution of equation 
© is: 



/ (y) = exp (a + by) 



(10) 




(12) 



which is numerically solved to obtain /o, or alternatively 
a. In reference 9], the solution of the equation Q with 
the boundary condition |(3J) is treated in detail. 

In figure^we plot the magnetization along the z-axis, 
S^, for the values of L = 3, = 2.5 • lO'^ (both in J 
units). This is obtained by inverting the above changes 
to the original angles. 

We calculate the total magnetization that appears in 
the d-spins system along the z-axis. The result is show 
in table m As can be seen, the existence of the magnetic 
distortion outside the localization region of the charge 
carrier leads to a partial compensation in the total mag- 
netization. 

With the analytical solution one can easily calculate 
the radius of the magnetic distortion. As can be seen 
from the solution, equations EJ, the radius is not 
properly defined, because the distortion created by the 
conduction electron in the magnetic system of d-spins 
reaches the complete chain. We choose to define the ra- 
dius as the distance at which the straight line with slope 



equal to 



de 

dy 



and passing throughout the origin at 



^0, reaches the value tt/2. This underestimates the ra- 
dius, but has the obvious advantage that in this way, the 
radius depends only on the value of the constant a, or 
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alternatively on the angle 



(13) 



Als o note that the value of (?o ~ and therefore R « 
TTy/J/SK. We obtain R = 6.18981 (in lattice constant 
units) for K ^ 2.5 ■ 10"^ (in J units). 

With the analytical solution it is also easy to calcu- 
late the energy of the ferron. The equation for the en- 
ergy, equation , contains also the energy for the system 
of d-spins even in the case with no conduction electron 
present. This energy has a value of £^o = ^{J + K)N. 
We define the energy of the magnetic polaron as -E — £^o- 
As the transformation that we perform to obtain the dif- 
ferential equation lO is a canonical transformation, we 
can use it again to calculate the energy. Taking this into 
account we transform the equation ^ and make again a 
power expansion in the lattice constant: 



E — -E'corc 2 



dylj 



1 fdev 




where E'coro = J (1 cos 26*0) - 2i cos 6*0 - t [IJ . There- 
fore the polaron energy is: 



-Epoi — E — Ef) — Ec, 



V8JK {l-sm0a) (15) 



as expected from the result for the ferromagnetic domain 
wall. Also note again that the value of Oq is small and 
therefore i?poi ~ —2{L — J) — t + VSJK, not only in the 
case of wide-band, but also in the case of double exchange 
semiconductor. 

The motivation for the energy calculation is to demon- 
strate that the ground state energy of a bound ferron 
with a magnetization compensating region is much lower 
than the ground state energy of a bound ferron without 
such a long range distortion. In table ^ we present the 
result for the ground state energy of the magnetic po- 
laron calculated from the equation IjlSf) . We also show 
for comparison the result in the case without compensat- 
ing region. This latter was calculated solving the set of 
equations |(2Jl numerically, and imposing = 7r/2. As 
expected the energy coming from equation (|15l) is lower, 
meaning that a true bound ferron is much more stable 
than bound ferrous considered previously in the litera- 
ture, as it was once again anticipated by Nagaev in 'rj. 

As can be seen from the equation H15|l . the main part 
of the energy of the magnetic polaron is concentrated in 
the core, that means at the sites g — —1,0, where the 
charge carrier is trapped. The energy in the compen- 
sating region, that means outside the sites g = —1,0, is 
very small and positive. This seem to be in contradic- 
tion with the previous discussion, in which the presence 
of a compensating region was presented as energetically 
favored. The physical explanation of this behavior is the 





Ground state energy 


6*0 




With comp. reg. 


-3.60101 


0.22011 


1.000195 


Without comp. reg. 


-2.68811 


0.63688 


1.60791S 



TABLE I: Ground state energies (in J units), canting angle 
in the core, 9o, and total magnetization along z-axis (in S 
units), Mz, for a ID magnetic polaron with, and without 
compensating region, using L = 3, K = 2.5 • 10~^ (both in J 
units) . 



following. The main energy scale in the problem that 
is coupled to the magnetic ordering is L, the s-d inter- 
action, in a wide-band semiconductor. It tends to put 
the d-spins in the core as parallel to the z-axis as pos- 
sible. But the d-spins in the core are connected to the 
rest of the chain through the d-d exchange term. The 
role of the compensating region is to isolate the d-spins 
in the core from the rest of the chain. This allows the 
d-spins in the core to be parallel to the z-axis with a high 
gain in s-d exchange energy (-Ecoro is high and negative) 
and a little loss due to the perturbation of the antifer- 
romagnetic ordering on the remaining part of the chain 
(the energy of the compensating region is positive but 
small). To better explain this point, we also show in ta- 
ble the canting angle of the d-spins of the core, for 
a bound ferron with compensation region, and a bound 
ferron without compensating region. As can be seen, the 
presence of the magnetization compensation region with 
the structure described above strongly reduces the value 
of canting angle. The role of the anisotropy is simply to 
set a scale for the radius of the magnetic distortion in the 
d-spin system. 

Apart for the solution treated here, there is another 
possible solution for the d-spin structure ,23]. This sec- 
ond solution is characterized by the absence of extended 
spin distortions. It corresponds to rotate by tt the d- 
spins of one half, say the negative one, g < 0, oi the 
chain. This make that the magnetization in the trapping 
region will be directed along x-axis, instead of z-axis. 
Mathematically, it corresponds to choose a different sym- 
metry condition for the d-spin system, namely — 0, 



and 5"^ 



5.. 



This solution saves more anisotropy 



energy and allows for a lower canting angle in the core, 
being at first instance the ground state in the range of pa- 
rameters chosen here. However, as the Heisenberg model 
has to be solved over a compactified ring, that is, the 
d-spin at g — > +00 is linked to the d-spin at g — > —00, 
one has to add an extra energy 2 J to the polaron energy, 
equation resulting that for the case of an isolated 

ferron the ground state energy always corresponds to the 
solution treated here. In the case of a finite density of fer- 
rous, this condition at infinity can be skipped and both 
solutions would be possible. The finite density case is left 
for further work. 

The three-dimensional case cannot be solved analyti- 
cally. We only give a simple estimate for the radius of the 
distortion for an antiferromagnetic coupling. We use the 
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three-dimensional analogue of the Hamiltonian . The 
case i^T' = is treated in reference 6J. To introduce the 
anisotropy we use a simple variational method. We as- 
sume that the distortion angle is given by the spherically- 
symmetrical asymptotic solution of Dc Gcnnes, 9 [r) ^ 
l/r^ for r < R, and that there is no distortion for r > R. 
The optimal value of R is obtained by minimization. For 
R large, we obtain R ^ {J / K) ^ . This provides an order 
of magnitude for the radius of the distortion. With this 
estimate for the radius of the distortion is now clear that, 
in real three-dimensional materials and over experimen- 
tal doping ranges, the (one-electron) ferrons treated here 
must overlap, forming magnetic droplets with a number 
of conduction electrons larger than one, as reported in 
the experiment. Therefore it is expected that the main 
conclusions of this article hold for the three-dimensional 
case, although further calculations are needed. 

To summarize, we have found the detailed structure 
of the one-dimensional d-spin system in the region sur- 
rounding of a bound ferron, completing the previous re- 



sults of 0, 0, • The main result is the appearance of 
a new length scale, namely the extent of the magnetic 
distortion created by the charge carrier. The existence of 
this distortion makes the ferron more stable. This may 
determine, together with the Coulomb interaction, the 
spatial distribution of magnetic droplets and their cou- 
pling in antiferroniagnetic semiconductors, such as un- 
derdoped manganites. Also it is related to the onset of 
the electronic phase separation at the very low doping 
range observed in these compounds. 
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